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ABSTRACT 

We present IRAC/MIPS Spitzer observations of intermediate- mass stars in 
the 5 Myr old A Orionis cluster. In a representative sample of stars earlier than 
F5 (29 stars), we find a population of 9 stars with a varying degree of moderate 
24/im excess comparable to those produced by debris disks in older stellar groups. 
As expected in debris disks systems, those stars do not exhibit emission lines in 
their optical spectra. We also include in our study the star HD 245185, a known 
Herbig Ae object which displays excesses in all Spitzer bands and shows emission 
lines in its spectrum. We compare the disk population in the A Orionis cluster 
with the disk census in other stellar groups studied using similar methods to 
detect and characterize their disks and spanning a range of ages from 3 Myr to 
10 Myr. We find that for stellar groups of 5 Myr or older the observed disk 
frequency in intermediate mass stars (with spectral types from late B to early 
F) is higher than in low mass stars (with spectral types K and M). This is in 
contradiction with the observed trend for primordial disks evolution, in which 
stars with higher stellar masses dissipate their primordial disks faster. At 3 Myr 
the observed disk frequency in intermediate mass stars is still lower than for low 
mass stars indicating that second generation dusty disks start to dominate the 
disk population at 5 Myr for intermediate mass stars. This result agrees with 
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recent models of evolution of solids in the region of the disk where icy objects 
form (> 30 AU), which suggest that at 5-10 Myr collisions start to produce large 
amount of dust during the transition from runaway to oligarchic growth (reaching 
sizes of ~500 km) and then dust production peaks at 10-30 Myr, when objects 
reach their maximum sizes (>1000 km). 

Subject headings: infrared: stars: formation — stars: pre-main sequence - 
open cluster and associations: individual ( A Orionis cluster) — protoplanetary 
systems: protoplanetary disk 



1. Introduction 



Planets are thought to form in circumstellar disks of dust and gas that result from the 
collapse of rotating molecular cloud cores during the star forming process. These primordial 
optically thick disks evolve due to viscous processes by which t he disk is accreting mat erial 
onto the star while expanding to conserve angular momentum (IHartmann et al. II1998I ). In 
addition, dust gra ins in the disk are exp ected to grow in size and settle toward the mid- 
plane of the disk (IWeidenschilringj 119971), where plane tesimals and planets can be formed 
( IBackman &: Paresce I Il993l ; iKenyon &: Bromley! 120091) . The dispersa l of primordial disks 



operate less efficiently for stars with lower mas ses (Mu zerolle et al 



2003 



Sicilia-Aguilar et al. 



20051 : lLada et al.l 120061 : iCarpenter et al.l 120061 : IHernandez et al.l l2007al : Kennedy &: Kenyon 
20091 ) . Particularly, for low mass stars (K5 or later), 90% of the stars lose their primordial 
disks by about 5-7 Myr (e.g. lHaisch et al.ll200ll ; IHernandez et al.l 120081 ). while for objects in 
the mass range of Herbig Ae/Be (HAeBe; types B, A or F early) stars, the corres ponding; 



time scale for primordial disk dissipation is less than 3 Myr (IHernandez et al.ll2005l ). 



The HAeBe stars, which have primordial optically thick disks, are the precursors of in- 
termediate mass stars with debris disks (like Vega or (3 Pic). Since the estimated lifetime for 
circumstellar dust grains due to radiation pressure, sublimation, and Poynting-Robertson 
drag is significantly smaller than the estimated ages for the stellar systems, the gas-poor 
disks observed in Vega type objects must be replenis hed from a reservoir, such as sublima 



tion of comets or collis ions between parent bodies (jChenl 120061 : IKenyon fc Bromley! 12004 



Dominik fc Decinl 120031 ) . Most obser vational stud i es of the evolution of d e bris disks have 



focused on the long time scales (e . g iDecin et al 1 120031 ; iRieke et al.l 120051 : ISu et al.l 12006 



Siegler et al.l 120071 : ICarpenter et al.l 120091 ). These studies show that the luminosities of 
the debris disks decay following a simple power law (with characteristic time-scales of sev- 
eral hundred millions yea rs), as the debris reservoi r of collisional parent bo dies diminishes 
and the dust is removed (IKenyon fc Bromley! 12004 : iDominik fc Decinl 120031 ). Observations 
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suggest a peak in the debris disk phenome non in intermediate mass stars at 10-30 Myr 



( [Hernandez et al.l 120061 ; ICurrie et al.l l2008al ); at the peak, debris disks are more frequent 



cal models of evolution of solids in the c 
disks may trace active planet formation 


isk f 


Kenvon & Bromlev 


2004, 


200d). Since debris 


Greenberg et al. 


1978: : 


3ackman & Paresce 


1993; 


Kenvon & Bromlev 


2008: 


Wvatt 


2008: 


Cieza 


2008 


), studies of debris disks younger than 



10-30 Myr are crucial to understanding the processes in which primordial disks evolve to 
planetary systems. 

The A Orionis cluster (also known as Collinder 69) represents a unique laboratory for 
studies of disk evolution because it is reasonably near and relatively populous, making possi- 
ble statistically significant studies of disk pro perties in a wide range of stellar masses. Based 



on main sequence fitting of early-type stars, iMurdin &: Penstonl ( 119771 ) obtained a distance 
of 400 ±40 pc and age of 4 Myr . Similarly, using Stromgren photometry and main sequence 
fitting. iDolan fc Mathie u (2001) derived a distance of 450±50 pc and a turnoff age of 6 Myr. 
Recently, iMayne fc Naylorl (120081 ) reported that the A Orio nis cluster and the a Orionis clus 



ter have an age of 3 My r. However, their disk frequencies (jBarrado y Navascues et al. 



2007 



Hernandez et al.l l2007al ) and comparison between empirical isochrones ([Hernandez et al. 
20081 ) suggest that the A Orionis cluster is older than the a Orionis cluster. For consis- 



tency with previous stud ies of disk and stellar population, we assume a distance of 450 pc 
and an age of 5 Myr (see iMathieul 120081 ) . 



The unprecedented sensitivity and spatial resolution provided by the Spitzer Space Tele- 
scope in the near- and mid-infrared windows are powerful tools to expand significantly our 
understanding of star and planet formation processes. In this contribution, we analyze the 
near- and mid-infrared properties of intermediate mass stars in the A Orionis cluster. This 
paper is organized as follows. In SJ2] we describe the observational data and the selection of 
intermediate mass star in the cluster. In §3] we present the method for identifying stars with 
infrared excesses. In SQl we compare the disk population of the A Orionis cluster with results 
from other stellar groups inferring the nature of disk emissions around intermediate mass 
stars. Finally, we give our conclusions in §51 



Observations 



2.1. Spitzer observations 



We have obtained near-infr ared and mid-infra red data of the A Orionis cluster using the 
InfraRed Array Camera(IRAC, iFazio et al.l 120041 ) with its four photometric channels (3.6, 
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4.5, 5. 8 fc 8.0 /xm) and t he 24/im band of the Multiband Imaging Spectrometer for Spitzer 



(MIPS iRieke et al.ll2004l ) on board the Spitzer Space Telescope. These data were collected in 
mapping mode during 2004 October 11 (IRAC) and March 15 (MIPS) as part of G. Fazio's 
Guaranteed Time Observations in Spitzer program 37. 

The IRAC observations were done using a standard raster map with 280" offsets, to 
provide maximum areal coverage while still allowing ~20" overlap between frames in order 
to facilitate accurate mosaicking of the data. Images were obtained in high dynamic range 
(HDR) mode, in which a short integration (1 s) is immediately followed by a long integra- 
tion (26.8 s). Standard Basic Calibrated Data (BCD) products from version S14.0.0 of the 
Spitzer Science Center's IRAC pipeline were used to m ake the final mosaics. Post-BCD data 



treatment was performed using custom IDL software (IGutermuth et al.ll2004l ) that includes 
modules for detection and correction of bright source artifacts, detection and removal of cos- 
mic ray hits, construction of the long and short exposure HDR mosaics, and the merger of 
those mosaics to yield the final science images with a scale of 1".22 pixel -1 . Since mosaics of 
channels 4.5 and 8.0 /im have a 6.5' displacement to the north from the mosaics of channels 
3.6 and 5.8 /im, the region with the complete IRAC data set (57.6'x54.6'; hereafter IRAC 
region) is smaller than the field of view (FOV) of individual mosaics (57.6'x61.2'). IRAC 
point source photometry was obtained using the apphot package in IRAF, with an aperture 
radius of 3". 7 and a background annulus from 3.7 to 8". 6. We adopted zero-point magnitudes 
for the standard aperture radius (12") and background annulus (12-22". 4) of 19.665, 18.928, 
16.847 and 17.391 in the [3.6], [4.5], [5.8] and [8.0] bands, respect ively. Aperture corr ections 



were made using the values described in IRAC Data Handbook (jReach et al. II2006I ). Final 



photometric errors include the uncertainties in the zero-point magnitudes (~0.02 mag). 

MIPS observations were made using medium scan mode with full-array cross-scan over- 
lap, resulting in a total effective exposure time per pointing of 40 seconds. MIPS observations 
cover more than 99% of the IRAC region. The images were processed using the MIPS in- 
strument team Data Analysis Tool (DAT), which calibrates the data and appli es a distortion 



corre ction to each individual exposure before combining it into a final mosaic (IGordon et al. 



20051 ). A second flat correction was also done to each image using a median of all the images 
in order to correct for dark latents and scattered light background gradients. The second 
flat correction is determined by creating a median image of all the data in a given scan leg, 
and then dividing each individual image in that scan leg by that median. Bright sources and 
extended regions are masked out of the data before creation of the median, so the nebulosity 
seen in the map has a negligible effect on any noise that might be added by the second flat. 
We obtained point source photometry at 24 /iin with IRAF/ daophot point spread function 
fitting, using an aperture size of about 5". 7 and an aperture correction factor of 1.73 derived 
from the STinyTim PSF model. The aperture size corresponds to the location of the first 
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airy dark ring, and was chosen to minimize effects of source cro wding and background in - 
homogeneities. The absolute flux calibration uncertainty is ~4% (jEngelbracht et al.l 120071 ). 
The photometric uncertainties are dominated by the background/photon noise. Our final 24 
/iin flux measurements are complete down to about 0.8mJy. 



2.2. Sample selection and optical spectra 



Our procedure for identifying the disk populati on around intermediate mass stars in 
the A Orionis cluster is similar to that described in iHernandez et al.l (120061 ). To find the 
bright stellar population of the A Orionis cluster, we selected from the Two Micron All Sky 
Survey (2M ASS) catalog^ stars with J<11. This limit was c alculated using the V magnitude 



( iCox 1120001 ) and the V-J color (IKenvon fc Hartmann 1995 ) for a main se quence star with 



spectral type F5 located at 4 50 pc (jPolan fc Mathieul 12001 



visual extinction of A v ~0.4 (IDiplas fc Savage 



Mathieul 120081 ) and assuming a 



199J). The reddening in the 2M ASS bands 



was calculated from the assumed visual extinction using the ICardelli et al.l (119891 ) extinction 
law with Rv—3.1. Our preliminary sample includes 159 stars located in the IRAC region. 

We found two relatively bright stars (V<12) without 2MASS counterpart when com- 
paring our preli minary sample with photometric studies of the b r ight population in the A 
Orionis cluster (bolan fc Mathieul l200ll : iMurdin fc Penstoni fl9771 ; buerr et all Il982h . One 
star (HP36862) is a close companion (<3") of the central bright star of the cluster (the A Ori- 
onis star), not resolved in the IRAC/MIPS images. The other star, HP245168, was included 



previously as member of the A Orionis cluster with spectral type B9 (IPolan fc Mathieull2001 



Murdin fc Penstoni 1 19771 ). Visual inspection of the 2MASS images reveals that HP 245168 
is a stellar source not included in the 2MASS catalog. We estimated 2MASS photometry 
for this source using differential photometry. We use 2MASS sources located inside 6' of 
radius to get the median value of differential photometry between the 2MASS magnitudes 
and the instrumental magnitudes obtained on the 2MASS images using an aperture radius of 
3" and a background annulus from 3" to 7". Our estimation for HP 245168 indicates 2MASS 
photometry of J=9.50, H=9.52 and K=9.4 5 with a pho t ometr ic error of 5%; we added this 
star to our preliminary sample. Using the iKharchenkol (1200 ll ) catalog (limiting magnitude 
of V=12-14) and visual inspection on the IRAC 3.6/zm image, we did not find additional 
bright stars not included in the 2MASS catalog. 

Figure [1] illustrates the procedure to select early type candidates (F5 or earlier) of the 
A Orionis cluster. The NIR color-color diagram (upper-left panel) for the preliminary sam- 



^izier Online Data Catalog, 2246 (R.M. Cutri et al. 2003) 
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ple r eveals that most of the stars are located on the main- sequence locus f lBessell Brett 



19881 ) indicating that, in general, these stars have little extinction and no NIR excesses. 



V1271 Ori), a well known Herbig Ae star with spectral type A5 (e.g., 


Wade et al.l 


20071; 


A eke et al. 


2005; 


Hernandez et al. 2004: de Winter et al. 2001 


; Herbst & Shevchenkol 


1999; 


Finkenzeller & Mundt 


1984 


), appears in the HAeBe star loci (I 


lernandez et al.ll2005t). where 



the NIR excess can be explained by emission from an optic ally thick primordial disk wit h 



a sharp dust-gas transition at the dust destruction radius (IDullemond fc Dominiki 12004] ) . 



The other star (HP 36 881) is a spectroscopic binary with a peculiar spectral type B9III 



( IPedoussaut et al.lll988l ). Since HD36881 does not lie near the ZAMS (eve n allowing for pho- 



tome tric errors or binarity), probably this star is a foreground nonmember (IDolan k Mathieu 



200ll . see also lower panel of Figured]). The dotted horizontal line in the NIR color-color di- 
agram is the [J-H] limit corresponding to a F5 star with Av=0.4. We identified 35 early type 
candidates located below this limit. We also included the Herbig Ae (HAe) star HD245185 
as a early type star of the A Orionis cluster. 

The upper-right pan el shows the vector point diagram for the early type candidates. 
Proper motions are from iRoser et al.l (120081 ). In this panel, we also show distributions of 
proper motions which can be represented by Gaussians centered at fi a * cos(5) ~ 0.0 mas 
yr _1 and fi$ ~ 0.3 mas yr -1 , with cr=2.5 mas yr -1 . Most of the stars are located in a well 
defined region represented for the 3a criteria circle in the vector point plot; however, the 
proper motion of the star A Orionis does not agree with the general trend of the cluster. 
To explain t hese peculiar proper motions (as other peculiar kinematic properties), s everal 
authors (e.g.. lMaddalena k Morrisll987l : ICunha k Smithlll996l : [Dolan k Mathieulll999[ ) have 
invoked the presence of a massive companion of the A Orionis star that became a supernova 
(SN), affecting the A Orionis's kinematics and removing nearby molecular gas at the ce nter 
of the cluster about 1 Myr ago when the supernova exploded (IDolan k Mathieull200ll ). If 
this scenario is correct, then the SN explosion may have affected the kinematics of other 
members of the cluster specially if those members were part of a multiple system, A Orionis 
+ SN, or were still in formation and embedded in their molecular cloud. So we can not use 
kinematic properties alone as a reliable criteria for membership. 

The lower panel in Figure [TJ shows the color magnitude diagram, V versus V-J, for 
the early type candidates. Mos t stars are locate d near or in the region defined by the 
ZAMS and the 5 Myr isochrone (ISiess et al.ll2000l ). Since it is well known that theoretical 
isochrones not includin g the birth l ine d o not accurately match empirical isochrones for in- 
termediate mass stars (jHartmannl 120031 ) . and the oretical evolutionary m odels often ignore 
some hydrodynamic processes in stellar models (IMamajek et al.l 120081 ). these theoretical 
isochrones are plotted as reference only and are not used for our selection of the photo- 
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metric candidates. However, comparing empirical isochrones in the low mass star range, 
we showed that the A Orionis cluster is in sim ilar evolutionary stage as other groups with 
ages quoted as 5 Myr ([Hernandez et all 120081), which is in agre ement with previous results 
(e.g., iBarrado v Navascues et al.l 120071 ; bolan k MathieulhoO^ . Two stars (HD36881 and 



HD37148) are located above the general trend described by the sample . The star HD36881 
was reported previously as a foreground star (IDolan fc Mathieull200ll ) as discussed above, 



and HD37148 has the largest proper motion of the sample (/x,5=-33.5 mas yr -1 ; see upper- 
right panel), even larger than the proper motion of the star A Orionis. We excluded those 
two stars from being members of the cluster. In this contribution, we also excluded stars 
located right-ward from the dotted line, which represents the limit corresponding to a F5 star 
with Av=0.4. Our final sample includes 29 early type candidates of the A Orionis cluster. 

Stars with detected infrared excesses (see §3D are encircled in the lower panel of Figure 
[TJ Optical spectra were obtained for these stars using the 1.3m McGraw-Hill telescope of 
the MDM Observatory equipped with the Boiler & Chivens CCD Spectrograph (CCDS). 
We used the 150 grooves per mm grating centered at 5300A. This configuration provides 
a nominal resolution of ~5A with spectral coverage from 3800A to 7100A. In Figure [2] we 
show spectra for the stars with infrared excesses sorted by spectral types. Spectral types were 
ob tained using the SPTCLASS tool I, an IRAF/IDL code based on the methods described 



m 



Hernandez et al.l (120041 ) . The HAe star HD 245185 sho ws emission lines sug gesting that 



magnetospheric accretion processes are present in this star (iMuzerolle et al.ll2004j ). The other 
stars do not show emission lines in their optical spectra indicating that magnetospheric 
accretion has stopped; this characteristic is typically observed in diskless stars and stars 
bearing debris disks . Mos t stars not observed with the CCDS have spectral types from 



Cannon fc Pickering! (119931 ) . For the stars without spectral type information, we estimated 



the latest spectral type they may have by assuming no reddening and interpolating the 
observed V-J color in the table of standard colors given by lKenyon fc Hartmannl (119951 ) . 



Table [T] shows optical and Spitzer data for the early type candidates. Columns (1) and 
(2) show the name and the 2MASS denomination of the stars. Column(3), (4), (5) and (6) 
show the IRAC magnitudes. Column (7) shows the 24 /im MIPS fluxes. The typical errors 
in the IRAC and MIPS magnitudes are 0.02-0.05 and 0.03-0.06, respectively. Column (9) 
shows the excess ratio of the measured flux at 24/im from that expected from the stellar 
photosphere (see §3]). Spectral types and references are in columns (10) and (11). Visual 
magnitudes and references are in columns (12) and (13). The last column shows the disk 
type around each star (see §3] and §1]). 



2 http:/ /www. astro. lsa.umich.cdu/^hernandj/SPTclass/sptclass. html 
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3.1. 



3. Disk diagnostic 
Identifying stars with infrared excess 



To characterize the disk population among stars of the A Orionis cluster, we need to 
identify photometric candidates that exhibit excess emission at the IRAC/MIPS bands. 
Since a disk produces greater excess emission above the stellar photosphere at longer wave- 
lengths, we selected the bands at 8/im and 24/im to identify stars with infrared emission 
due to the presence of disks. First, we identified stars with 24/im infrared excess above 
the photospheric level in the upper panel of Figure [3j In the upper portion of this panel, 
we plot the K-[24] color distribution for the sample. The distribution of stars around the 
expected photospheric color (K-[24]~0) describes a Gaussian centered at K-[24]=-0.09 with 
cr=0.16. The 3cr boundaries (vertical dotted lines) represent the photospheric colors; thus 
stars with excess at 24/im have colors K-[24]>0.4. Besides the HAe star, which displays 
the largest excess at 24/im in our sample, we detected 9 stars (hereafter debris disk can- 



in other young stellar re gions (e.g 



Youne; et al. 


2004; 


Hernandez et al. 


2006; 


Gorlova et al. 



20071 ; ICurrie et al.ll2008d .lbl). and with no emission lines in their optical spectra (see Figure [2]). 
As referenc e, we displayed the loc a tion o f the spectral type sequence, using the standard V- J 
colors from lKenyon fc Hartmannl (119951 ). and the corresponding K-[24] photospheric colors 
from the STAR-PET tool of Spitzer Science Center^ 

The lower panel of Figure [3] shows the [3.6]-[8.0] color - [3.6] magnitude diagram illus- 
trating the procedure to detect infrared excess at 8/im. The photospheric limits are defined 
using a 3a criteria based on the [3.6]-[8.0] color distribution, which can be represented by 
a Gaussian centered at [3.6]-[8.0]~0 with o=0.06. The error bars include the saturation 
effect, that were estimated using the deviation from the expected photospheric colors for the 
brightest objects detected in each IRAC band. Thus, the apparent 8/im excess observed in 
stars with [3.6] <6 mag are produced by saturation effect in the [3.6] band. Besides the HAe 
star which displays excesses in all IRAC bands, the star HD 245370 shows small excess at 
8/im; this star also shows the third largest excess at 24/im ( K-[24]~ 2). 



3 http: / /ssc. spitzer. caltech.edu/tools/starpet/ 
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3.2. Spectral Energy Distributions 



We have plotted the SEDs for the disk bearing stars found in our early type sample of 
the A Orionis cluster. Figure H] shows the SEDs for the Herbig Ae star HD 245185 and for 
the debris disk candidates sorted by spectral types. The color excess, E V -j, was calculated 
using the observed V- J color and the st andard color obtained i nterp olating the spectral type 
(see §2.2p in the calibration given by iKenyon fc Hartmannl (119951). Using the Ey_j and 
the extinction relation for normal interstellar reddening ( jCardelli et al.lll989l ). we obtained 
the visual extinctions reported in Figure HI Since the J magnitude of the HAe star could 
have contribution f rom the disk, we calculate the Ay for this star using the B-V color from 
Kharchenkd (120011 ) . Out of 10 stars plotted in Figure HI two stars (20%) have Ay larger 
than the value used as cut-off in §2.21 (Av=0.4). Since stars bearing disks are expected to 
have larger extinctions produced by the circumstellar material, we can estimate that the 
completeness in our selection of stars earlier than F5 is larger than 80%. A F0 star needs to 
have A v >0.8 to be rejected from our sample (see Figured]). 

Dotted lines in these plots represent the s tandard colors (or photospheric levels) for a 
given spectral type ( Kenyon fc HartmamJl995 ). Figure H] also shows the excess ratio at 24/^m 
( IRieke et all 120051 ) . which was calculated using the relation: E 2 4 = lo - 4 *^ - ! 24 !" 1 " ' 09 ), where 
-0.09 is the median value of K-[24] for the photometric candidates with MIPS detections (see 
§3.ip . In summary, we detected 10 early type disk-bearing stars in the A Orionis cluster. The 
HAe star HD 245185 is the only star that shows infrared excesses in all the Spitzer bands, in 
agreement with the presence of a primordial optically thick disk; emission lines in its optical 
spectrum indicates that HD 245185 is accreting material from the disk to the star. Only the 
debris disk candidate HD 245370 shows infrared excess at 8/xm, the remaining debris disk 
candidates only exhibit infrared excess at 24/mi. 



4. Debris disk frequency in young stellar populations 



In iHernandez et al.l (120061 ). we reported that debris disks in the range of ages from 2.5 
to 150 Myr are more frequent and have larger 24/nn excess around 10 Myr. This peak in the 
debris disk phenomenon is in agreement with models of evolution of solid s in the region of the 
protoplanetary di sk where icy planets are predicted to form (30-150 AU; IKenyon fc Bromley 
2008T ). Obs ervational results from the double cl uster h and y Pers ei (10-15 



2004 



2005 



Myr; ICurrie et al.ll2008af ) and the cluster NGC2232 (25 Myr; ICurrie et al.l 12008 bh support 
our finding. However, it is still an open question whether the disk emission observed at 
24/zm in debris disks candidates younger than 10 Myr originates from remaining primor- 
dial dust or from second generation dust produced by collisional cascades in the disk (e.g. 
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Hernandez et all 12009 : Hillenbrand! 120081 : IWvattll2008l : ICiezal 120081 ). To tackle this question, 
we have compared the disk population in the A Orionis cluster to disk populations of several 
stellar gr oups located within 500 pc of t he Sun detected and characterized with the same 
methods JHernandez et alJbood . [2007al lbl. 1 20081 . Hernandez et al, in preparation). 



Figure [5] shows the theoretic al stellar fluxes a t 24/im for diskless stars in the stellar 
populations of Table [2j Using the ISiess et al.1 (120001 ) evolutionary tracks and the properties 
listed in Table [21 we estimated the apparent K magnitude for stars in each stellar group. The 
24 jum fluxes were calculated using the corresponding K-[24] photospheric colors obtained 
from the STAR-PET tool of Spitzer Space Center and the zero point magnitude defined 
for MIPS observations (~ 7.17 Jy@). Given the limiting flux at 24 /iin (0.5-0.7 mJy), and 



using the ISiess et al.l (120001 ) evolutionary tracks to estimate the K magnitude for stars in 
each stellar group, we inferred that we can detect any excess at 24 /im in stars F0 or earlier 
with low or moderate reddening. Figure [5] shows the limit for diskless stars (E24=l). The 
detection thresholds at 24//m for intermediate mass stars in different clusters are given by the 
width of the photospheric region defined by the K-[24] color distributions. These detection 
thresholds are nearly the same in all stellar populations of Table El in which we can detect 
24/zm excess in intermediate mass stars with E 2 4 > 1.6. We also show in Figure the 
detection threshold for stars bearing disks with E 2 4=5. Figure [5] indicates that for all our 
surveys the spectral type limit for detecting stars bearing disks with E24 >5 is Ml±l. 

Using the same method for 8 /im, we estimated that a diskless star in our surveys with 
spectral type M5 has [8.0] <13. This is well above of the IRAC detection limits of our surveys 
([8.0]~14.0-14.5). The detection threshold at 8/im for a M5 star is defined by the error of the 
IRAC SED slope. The IRAC SED slope, used to identify infrared excess at 8/im in low mass 
stars, is calculated from the [3.6]- [8.0] color and its error is dominated by the photometric 
error at [8.0]. Thus, using the typical error in the IRAC SED slope at [8.0] =13 (~ 0.25 for 
all surveys), we can detect disks around stars with spectral type M5 or earlier with a excess 
ratio at 8 /im larger than 1.2. 

Figure [6] shows observed disk frequencies as a function of color/spectral type for several 
young stellar populations, the clusters o Ori, A Orionis, 7 Velorum, a region in the Orion 
OBlb subassociation, and the stellar aggregate 25 Ori. Frequencies were calculated for 
several bins of spectral types in each stellar population. 

The bin corresponding to the intermediate mass stars includes stars ranging from 
B8 to F0; most of the intermediat e mass stars with infrared excesses are in this range 
( jHernandez et al.ll2006l . l2007al 120081 ) . In general, the intermediate mass stars of the stellar 



4 http:/ /ssc. spitzer. caltech.edu/mips/calib/ 
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populations plotted in Figure [6] have spectral type information available in the literature 



(iHouk 



1978 



Hernandez et al 



2005 



Houk & Swift 1999; Nesterov et al. 11995: Kharchenko 2001; Caballero 2007 



20(3); optical photometric colors were used to estimate spectral types 



for stars without spectral type information. In agreement with the extinctions listed in Table 
121 the reddening for intermediate mass stars in these stellar groups are relatively low. Except 
for the star HD36219 located in the Orion OBlb subassociation (with Av=1.46), all stars in 
the intermediate mass samples plotted in Figure [6] have Ay<l- 

Disks around low mass stars were split in three bins of spectral types: from K0.5 to M0. 5, 
from M0. 5 to M3.5 and from M3.5 to M6 (hereafter; the low mass bins). Given the low red- 
dening expected for the stellar populations plotted in Figure O photometric colors represent 
a good approximation of spectral types when spectroscopic informati on is not available for 
the st ellar groups. Thus, we used the standard table of colors given by lKenyon fc Hartmann 
(119951 ) to obtain spectral type estimations. Observed disk frequencies of low mass stars in 
the A Orionis cluster were calc ulated analyzing the IRA C and MIPS photometry of member s 
confirmed spe c trosco pically by iDolan &: Mathieul (120011 ) ; iBarrado y Navascues et al.l (120071 ) ; 



Maxted et al.l (120081 ) . and ISacco et al.l (120081 ) using radial velocity distributions and the 
presence of Li I in absorption (Hernandez et al in preparation). Spitzer point source pho- 
tometry for these low mass members was obtained using the same methods d escrib ed in 
£ 12.11 Disk detection follow the method described in Hernandez et al. rl2007al lbl. [2008h . We 
used NI R photometry from the 2 MAS S catalog (jCutri et al.l 120031) and optical photome- 



try from IDolan fc Mathieul (120021 ) and IBarrado y Navascues et al. 



(12007( 1 . Since the R-J 



color is more complete for the spectroscopic members compiled for the A Orionis cluster, we 
used this color to calculate photometric spectral types and separate the sample in the bins 
described above. For comparison, we plot ted the disk frequencies obtained from the disk 



census m 



Barrado v Navascues et al.l (120071 ). Since the IRAC photospheric region defined in 



Barrado y Navascues et al.l (120071 ) does not take into account photometric errors, the num- 
ber of stars bearing disks could be overestimated in the lowest mass bins ( Maria Morales 
Calderon, private communication). 



For the a Orionis cluster, we used the disk census in IHernandez et al.l (j2007al ). To 
reduce the contamination by non-members, observed dis k frequencies of low mass stars were 
calculated usi ng the spect r oscopi c me mbers compiled in IHernandez et al.l (j2007al ) , members 



confirmed by Sacco et al. ( 2008 ) and Maxted et al. ( 20081 ). and the X ray members from 



Franciosini et al.l (120061). Since most of th is sample does not have sp ectral types, we used 
the c olor V-J in IHernandez et al.l (j2007al ) and the standard colors (IKenyon fc Hartmann 
19951 ) to separate the low mass sample in the low mass bins. 



For the Orion OBlb subassociation and the 25 Orionis aggregate we used the disk census 
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in IHernandez et al.l (l2007bl ). Observed disk frequencies were calcula ted using confirmed 



mem bers and spectral types calculated using our SPTCLASS code (jBricefio et al 



2005 



20071 . Briceno et al in preparation). Since our samples of confirmed members in the Orion 
OBlb subassociation and the 25 Orionis do not go as early as K0.5, the first low mass 
star bin in these groups does not cover the same spectral type range as the other stellar 
populations in Figure [6] (from K2 for the 25 Orionis aggregate and from K4 for the Orion 
OBlb subassociation). 



Finally, for the 7 Velorum cluster we used the disk census in IHernandez et al.l (120081 ) . 
Observed disk fre quencies were calculate d using the number of members expected as a func- 
tion of V-J color (IHernandez et al.ll2008l ); this color was used to separate the sample in the 
low mass bins. Using empirical Li-depletion studies and empirical isochrone comparisons, 
the 7 Velorum clu ster appears to be slightly older than the A Orionis cluster (5-7 Myr; 



Jeffries et al.l 120091 ) and younger than the 25 Orionis stellar aggregate (IHernandez et al. 
2008h. 



In Figure El the observed disk frequencies in the low mass stars bins were calculated 
using the infrared excess at 8 /im. Since a disk produces greater emission above the stellar 
photosphere at longer wavelengths, those stars with excess at 8 jum are expected to have 
excess at 24 /im. However, since the sensitivity of the 24 /im MIPS band is lower than that 
of the IRAC 8 /im channel, some of the faintest stars bearing disks are below the MIPS 
detection limit. We also included as stars bearing disks, stars with small or no infrared 
excesses in the IRAC bands but detectable excess at 24 /im (transitional disk candidates). 
The transitional disk candidate s represent a relative small fraction of the disk - bearing stars 



in th ese stellar groups (<10%; IHernandez et al.l l2007bl ; lErcolano et al.l 120091 ; lUzpen et al. 
20091 ) . so they do not significantly affect our conclusions. 



The infrared excesses observed in low mass stars are unlikely from second generation 
dust. Spitzer observations of 314 solar- type starj§ in t he Formation and Evolution of Plane- 



tary Systems (FEPS) Legacy program (jCarpenter et al.ll2009f ) indicate that the dust temper- 
ature of debris disks (Td ust ) ranges from ~50 K to ~200 K, with a fractional dust luminosity 
i f du.«t. = Lduxt./L *) less than 10~ 3 . Adopting these values, following the two-para meter model 
of lWyattl (120081 ). and using the stellar parameters for a K0 star at 3, 5 and 10 Myr (jSiess et al. 



200ol ). we estimated that debris disks around low mass stars show small or no excess at 8/im 
(Eg < 1.2) and modest excess at 24 /tm (E24 < 4). The two-parameter model includes two 
fundamental observable parameters of debris disks (T dust and fdust) and it is useful to es- 
timate the expected infrared excesses produced by a second generation disk with a unique 



5 ranging in spectral types from late F to early K and in the range of ages from 3 Myr to 3 Gyr 



13 



Td U st', however, notice that we need a more detailed model to reproduce the SED that for 
some disks can only be explained by the presence of several dust populations with different 
temperatures. In general, low mass stars bearing disks in our census show excesses much 
larger than those expected for debris disks (E 24 >25), thus the low mass bins in Figure M 
represent frequencies of primordial disks with little or no contamination by second generation 
dust. 



In contrast, the two-parameter model (jWyattl 120081 ) applied to debris disks around AO 
stars with dust temperatures ranging from 50 K to 200 K and fractional luminosity of 10~ 3 
indicates maximum E24 values of 29, 21 and 21 at ages of 3, 5, and 10 Myr, respectively. In 
general, debris disks around intermedia t e mass stars exhibit les s infrared excess than these 
maximum levels (e.g. iRieke et al.l 120051 ; iHernandez et al.l 120061 ; ICurrie et al.l l2008al ; IWyatt 



20081 ). Study ing the MIPS o bservations of 160 A- type stars (with age range from 5 Myr 
to 850 Myr), ISu et al.l (120061 ) show that debris disk systems have fd ust <10~ 3 and most of 
them (>90%)are characterized by dust temperatures between 50 K and 200K (although 
for stars older than 400 Myr the dust temperatures of the disks appear in a narrow range 
between 50 and 150 K). Since our debris disks candidates have similar levels of 24/xm excess 



than those reported in older debris di sks systems (e.g., IRieke et al.l 120051 ; ISu et al.l 12006 



Gorlova et al.ll2007l ; ICurrie et al.ll2008al ). we expect a similar rang e of T du.it and f^,,,^ D ebris 
d isks with Td„..<tt > 200 K (and /or fdust > 10~ 3 ) are very scarce ( Zuckcr man &: Song 



Willia ms fc Andrewsll2006l : ISu et al.ll2006l : Brvden et al. 2006: iRhee et al. 



2007 



Wvatt 



2004 



20081 : 



Uzpen et al.ll2009i ). This kind of massive debris disks (like HR4796A and /3 Pic) can be above 
our estimations of E 2 4 and could be explained if the observed dust is produced by eventual 



collis ions between two large objects (> 1000 km: iKenyon fc Bromleyil2004l . 120051 : ISong et al. 



20051 ). However, an alternative expla nation is that tho se objects are in a phase between 
HAeBe star and true debris systems (IRieke et al.l 120051 ) retaining some of the primordial 
dust. The star HD245370 shows excess at 8/xm that can be explained by the two-parameter 
model assuming a T dus t ~ 250 K, however this star is outside of the intermediate mass bin 
plotted in Figure (6) In our samples of intermediate mass stars, only three stars bearing disks 
show excess at 24/rni larger than the theoretical limits calculated above: HD 290543 in the 
Orion OBlb subassociation (i? 24 =537.5), HD 245185 in the A Orionis cluster (£ l 2 4=1129.3), 
and V346 Ori in the 25 Orionis aggregate (_E 2 4=376.4). These three stars also show Ha in 
emission indicating the presence of an accreting primordial disk. The remaining intermediate 
mass stars with disks exhibit characteristics normally present in debris disk systems, they 
have E24 <12.5 and no emission lines in their optical spectra. 

In summary, in the low mass bins the detection threshold is such that only primordial 
disks can be detected. In contrast, both primordial disks and debris disks can be detected 
in the intermediate mass bin. We may take then the observed frequency of disks in the 
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low mass bins in Figure E] as indicator of the frequency of primordial disks. Figure E] shows 
that the observed disk frequencies in K and M stars decrease as the stellar groups become 
older. Previous studies have shown that in young stellar populations pri mordial disks dis- 



sipate faster with high er stellar masses, for instance, in I C348 (2-3 Myr; lLada et al.l 12006 



Currie & Kenvonl 120091 ) . NGC 2264 f~3 ISung et all 120091 ). and a Ori in Figure El Accord- 
ing to this trend, the frequency of primordial disks around intermediate mass stars should 
decrease with age in proportion to the observed decrease in the primordial disk frequency 
in the low mass bins. However, Figure El shows that the observed frequency of disks in the 
intermediate mass bin increases after ~ 5 Myr. Since in this bin both primordial and debris 
disks can be detected, and since the primordial disk frequency is expected to decrease, the 
increase in disk frequency points to an increase in the frequency of debris disks. Assuming 
a binomial distribution for the errors of the disk frequencies, the significance of this increase 
ranges from 68% ( la) to 95%(2.0 a ) . Th e inc rease in frequency of debris disks has been 
reported before by IHernandez et al.l (120061 ) and ICurrie et al.l (j2008bl ) at ~ 10-15 Myr. The 
reversal in the observed disk frequency in the intermediate mass bin shown in Figure [6] shows 
that this increase is already present at much earlier ages, at ~ 5 Myr. Enough secondary 
dust must have been created by 5 Myr to make second generation dust detectable. This sec- 
ond generation dust must have been replenished from a reservoir, such as collision s between 



parent bodies (e.g.. IChenll2006l ; iDominik &: Decinll2003l ; iKenyon fc Bromley! 120041 ) 



Kenyon fc Bromley! (120081 ) described a suite of numerical calculations of planets growing 
from ensembles of icy planetesimals at 30-150 AU in disks around 1-3 M stars. These mod- 
els predict that the 24/zm excess from second generation dust rises at stellar ages of 5-10 Myr 
(during the transition from runaway to oligarchic growth, reaching sizes of ~500 km), then 
peaks at 10-30 Myr (when the first objects reach their maximum sizes, >1000 km) and finally 
slowly declines (as the reservoir of small sol ids in the disk diminishe s) . The expected theoret 



ical tr end has been observed previously by IHernandez et al.l (120061 ) and ICurrie et al.l (j2008a 



2009aJ). However, Figure E] suggests the second generation nature of the excesses observed at 



disks ( Kenyon &: Bromley 



2004 



2005 



to occur faster in the inner disk than the outer disks for primordial disks ( 


Weidenschilline 


1997; 


Dahm & Hillenbranc 


2007 


; lLada et al. 


2006; 


Sicilia-Aeuilar et al. 


2006 


) and for debris 



200 



J, there is a possibility that the debris disk can- 



didates in our studies have an outer primordial disk component which does not contribute 
significantly to the emission observed at A < 24 /im. An example of this kind of object is 
AU Mic, a ~ 12 old M-type star with a collisional evolved debris region in the inner part of 
the disk, and with ou ter regions composed by pristine material that is still part of the rem- 
nant primordial disk (IMetchev et al.ll2005l ). Additional observations at longer wavelength or 
observation of gas component of the disk are needed to confirm the complete dissipation of 
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the ou ter primordial disks in the young debris disks candidates. Finally, ICurrie fc Kenyon 
( 120091 ) suggest that debris disk s can be present in the 2-3 My r IC 348 cluster. Using theo- 
retical models of dust removal (ITakeuchi fc Artymowiczll200ll ). they estimate that the dust 
removal timescale is less than l/20th of the median age for IC 348 stars, even if there is some 
residual gas in an optically thin disk. Thus, we can not reject the possibility that some of the 
disks detected around intermediate mass stars of the a Orionis cluster are second generation 
disks. We need additional studies to confirm the actual nature of the intermediate mass 
stars bearing disks detected in this cluster. 



5. Summary and Conclusions 



We have used the IRAC and MIPS instruments on board the Spitzer Space Telescope 
to study the frequencies and properties of disks around the intermediate mass population 
of the 5 Myr old A Orionis cluster. We have selected 29 intermediate mass stars (~F5 
or earlier) using optical and 2MASS photometry. One star, HD 245185, displays excesses 
in all IRAC band, strong excess at 24 /im and emission line in its optical spectra, sup- 



dial disk ( 


Finkenzeller & Mundt 


1984; 


Hernandez et al. 


2004; 


Acke et al 


2005 



Wade et al. 



1984; Herbst fc Shevchenko 1999; de Winter et al. 2001 



20071 ). Additionally, we have identified 



nine stars as debris disks, which exhibit modest excess at 24/im and very small or no excesses 
in the IRAC bands. The SEDs of our debris disks sample are comparable to those observed in 
debris disk populations of older stellar groups. As expected in stars with second generation 
disks, these stars do not show detectable accretion indicators in their optical spectra. 

We have combined observations of the A Orionis cluster with those of other stellar groups 
of ages raging from 3 to 10 Myr located at < 500 pc from the sun, studied using similar 
methods to detect and characterize their disks. Because of observational thresholds, the 
observed disk frequencies observed in K and M stars reflect the primordial disk populations, 
while the observed disk frequencies in intermediate mass stars encompass both primordial and 
debris disks. As expected in primordial disk evolution the observed frequencies in K and M 
stars are smaller in older stellar groups. However, observed disk frequencies in intermediate 
mass stars, represented by a bin in spectral type from B8 to F0, show a behavior in sharp 
contrast with that expected for primordial disk evolution. In particular, the disk frequency 
in intermediate mass star increases from ~20% at ~3 Myr, to ~40% at 5 Myr and to ~50% 
at ~10 Myr. Since the timescale for primordial disk dissipation is dependent on stellar 
mass, the expected frequency of primordial disks around intermediate mass stars at ages 
> 5 Myr, as inferred from the corresponding frequency in low mass stars, is very low. In 
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stellar population of 5 Myr or older, the observed disk frequencies in intermediate mass stars 
are larger than in low mass stars. Thus, the increasing in the observed disk frequency at 
these ages indicates that these disks are dominated by second generation dust. This result 
agrees with models of evolution of solids in the disks (IKenyon fe Bromleyl 120081 ) . in which 
at 5-10 Myr collisions start to produce copious amount of dust during the transition from 
runaway to oligarchic growth. When oligarchs reach sizes ~500 km, collisions between 1-10 
km objects produce debris instead of mergers and the leftover planetesimals are ground to 
dust. 



We thank Thayne Currie and Maria Morales Calderon for insightful communications, 
and the anonymous referee for input that greatly improved this manuscript. This work is 
based on observations made with the Spitzer Space Telescope (GO-1 0037), which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology under a contract with 
NASA. J.Hernandez and N. Calvet gratefully acknowledge support from the NASA Origins 
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References: (1) ICannon fc Pickerind dl993h : (2) This work; (3) photometric spectral types; (41 iKharchenkol (|200lh ; (5) bolan fc Mathieul l|2002h 
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Table 2. Properties of stellar populations 
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*The last column represents the later spectral type in which we can detect infrared 
excess at 24/^m with E24 >5 (see Figure [5} 



Referenc es: ( 11 iBrown et al.l l ll998h; (21 ISherrv et al.l ||200S|); (31 iDolan fc Mathieul 

i200ll); (41 iMathieul (12008^; (5) iBriceno et alj i2007f); (6) Ivan Leeuwenl ll2007h ;C7 



iMillour et alj ll2007ll; (81 iNorth et al.l i2Q0j); (9) IBrown et al.Nl994h: flOl (Beiar et a.1 
l ll999h : (lll lDiplas fc Savagj lll994lH l2l IPozzo et al.l (1200011 7(131 1 Jeffries et alj J2009T) 
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Fig. 1. — Diagrams illustrating the selection of photometric candidates of the A Orionis 
cluster. The HAeBe star HD 245185 (triangle), the star Lambda Orionis (square), the stars 
HD 36881 (diamond) and HD 37148 (inverse triangle) are marked. The upper-left pane l 
shows the 2MASS color-color diagram. The standa rd sequences from Bessell &: Brettl (119881 ) 



are shown in dashed lin es. The locus of the CTTS (IMeyer et al.lll997l ) and the HAeBe stars 



( [Hernandez et al.ll2005l ) are defined by a long-dashed line and by a solid box, respectively. 
The photometric limit for a F5 star with reddening of Ay=0.4 is represented by a dotted line 
in the upper-left panel and in the bottom panel. We selected stars below the photometric 
limit. The upper-right panel shows the vector point diagram for stars selected in the upper- 
left panel. In general, these stars have similar proper motions; however, since the supernovae 
explosion could affect the kinematic properties of the cluster, kinematic criteria can not be 
used alone to reject probable member of the cluster. Stars with discrepant proper motion 
are plotted with crosses. Finally, the bottom panel shows the color-magnitude diagram V-J 
versus V. Our sample of early type star in the A Orionis cluster are stars located left-ward 
(bottom panel) the photometric limit. We rejected as member the stars HD 36881 and HD 
37148 which are located above the expected sequence. Large open circles represent stars 
with infrared excesses (see $3j) 
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Fig. 2. — Optical spectra of the intermediate mass stars bearing disks sorted by spectral 
types. The spectra were obtained with the CCDS spectrograph on the 1.3m Telescope of the 
MDM observatory. Spectral types were obtained using the SPTCLASS code, some features 
used for spectral classification are marked in the panels. Emission lines are detected only in 
the HAe star HD 245185 
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[3.6]-[8.0] 

Fig. 3. — V-J vs. K-[24] color-color diagram for intermediate mass stars in the A Orionis 
cluster (Upper Panel). Dotted lines define the locus expected for diskless stars based on the 
K-[24] color distribution. Nine stars exhibit moderate excess at 24/im (K-[24]<3) while the 
HAe star (HD 245185; triangle) exhibits strong excess at 24/xm (K-[24]~7.5). The bottom 
panel shows the color-magnitude diagram, [3.6] vs. [3.6]-[8.0], used to detect stars with excess 
at 8/xm. Dotted lines define the photospheric limit based on the [3-6]- [8.0] color distribution. 
The two brightest stars in this plot are affected by saturation at 3.6/im. Since these stars 
dot not show excess at 24/im, the excess detected at 8/xm is not real. Thus, in addition to 
the HAe star (triangle), we have detected one star with small excess at 8/im (HD 245370). 
The square represents the central star A Orionis 
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Fig. 4. — SEDs for intermediate mass stars with infrared excess. Each pane l show s the 
spectral type of the object (see §2.21) . and the excess ratio, E24 (IRieke et al.ll2005l ). cal- 
culated from the K-[24| co lor. Dotted lines show the corresponding photospheric colors 
( jKenyon fc Hartmannlll995l ). Except for the star HD 245185, Ay's were calculated using the 
V-J color. Since the J magnitude of HD 245185 can be contaminated by disk emission, we 
have used the color B-V to calculate the reddening in this HAe star. All SEDs are corrected 
by reddening and normalized at the J-band. 
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Fig. 5. — Theoretical stellar fluxes at 24/im. Combining the ISiess et al.l (120001 ) evolutionary 
tracks and the corresponding K-[24] photospheric colors from STAR-PET tool of Spitzer 
Science Center we estimated the theoretical flux expected for the stellar populations with 
distances, ages and extinctions listed in Table [2J Given the limiting flux at 24 /zm (0.6-0.7 
mJy), we plotted the limit where we can detect any excess at 24/im (E 2 4 =1). For reference, 
we also plotted the limit where we can detect stars bearing disks with E24 > 5. Vertical thin 
lines define the spectral type range in which we can detect stars bearing disks with E 2 4 > 5. 
The Arrow represents reddening vector of Ay=20. 
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Fig. 6. — Disk frequency versus spectral types for several star forming regions. In the 
stellar population of 3 Myr, disk frequencies for A-type stars are lower than disk frequencies 
around low mass stars in agreement with the expected trend for primordial disk dissipation. 
In contrast, at 5 Myr and older, disks around A-type stars are more frequent than in low 
mass stars indicating than second generation dusty disk must dominated the disk population 
at higher stellar masses. The low mass bins are centered at the spectral types K6.5, M2.5 
and M4.5. Using these values we calculated a reference detection threshold for each bin (see 
Figure EJ): -E24 ~ 3, 7, 10 (from left to right). The detection threshold for the intermediate 
mass bin are given by the width of the photospheric region in each stellar group (e.g. Figure 
E}, E 2 i ~1.6. For the A Orionis cluster, we plotted with X's the disk frequencies calculated 



fr om the disk census inlBarrado v Navascues et al.l (120071 ) . Since the IRAC thresholds defined 



